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Abstract 

We consider a quantum model of a nanomechanical flexing beam resonator interacting with a 
bath comprising a few damped tunneling two level systems (TLS's). In contrast with a resonator 
interacting bilinearly with an ohmic free oscillator bath (modeling clamping loss, for example), 
the mechanical resonator damping is amplitude dependent, while the decoherence of quantum 
superpositions of mechanical position states depends only weakly on their spatial separation. 



PACS numbers: 85.85.+j,03.65.Yz 



I. INTRODUCTION 



The past few years have seen dramatic progress towards achieving the necessary con- 
ditions for demonstrating macroscopic quantum behavior in mechanical systems.^ Schemes 
under investigation typically consider either micronscale mechanical resonators that are elec- 
trostatically coupled to superconducting qubits (quantum electromechanical systems) 2 * 3 * 4 * 5 ^ 
or larger mechanical mirror resonators that couple via radiation pressure to light trapped 
in an optical cavity (optomechanical systems).^ One of the prime motivations for demon- 
strating quantum behavior is to deepen our understanding of the so-called quantum-classical 
divide, in particular how classical dynamics emerges from the underlying quantum dynam- 
ics as system sizes (i.e., mass/energy content) increase P It is commonly accepted that the 
environmental degrees of freedom with which the mechanical resonator's mode of interest 
interacts is responsible for the emergence of classicalityJ^ In particular, the environment 
is thought to cause the rapid decoherence of initial quantum superposition states of the 
mechanical mode, resulting in an apparent classical mixture of the states. There is a con- 
siderable body of theoretical work investigating the effective quantum dynamics of open, 
single particle systemsP^ However, largely for reasons of calculational convenience, much 
of the effort has been devoted to the solvable model of an environment comprising non- 
interacting oscillators that are bilinearly coupled to a single oscillator system.^^ In light 
of the experimental progress mentioned above, an important issue is the actual nature of 
the dominant mechanical resonator mode environments. At the very low (i.e., cryogenic) 
temperatures to which the resonators must be cooled in order to observe quantum effects, 
it is not a priori obvious that the resonator mode dynamics can be mapped onto that of 
the oscillator system-oscillator bath model. In the present paper, we focus on a type of 
environment degree of freedom that is known to be relevant at low temperatures, namely 
the tunneling two level system .H 3 -!^ 

Tunneling two level system (TLS) defects were first invoked in the early seventies in or- 
der to account for the observed thermodynamic properties of amorphous, dielectric materials 
at low temperatures.^^ Further, convincing evidence for their presence was provided by 
acoustic phonon pulse decay and phonon pulse echo experiments) 17 * 18 ^ In particular, these 
experiments verified the characteristic saturation of TLS's with increasing acoustic pulse 
intensity for resonant phonon absorption and also measured TLS relaxation and dephasing 



2 



times. TLS's have recently received renewed interest as one of the main decay /decoherence 
mechanisms for superconducting qubitsP ^ 21 * 22 * 23 * 24 * 25 * 26 * 27 ! Signatures include inducing reso- 
nant splittings in the qubit energy level spectra! 2 - * 21 * 26 ^ and saturation of microwave power 
absorption by the dielectric oxide layer of the qubit tunnel junctions J 21 * 25 ! In contrast with 
the bulk amorphous dielectric materials involved earlier investigations, the much smaller, 
micronscale sizes of the superconducting qubits with gigahertz frequency energy level sep- 
arations exceeding the dilution fridge fcgT thermal energies point to a distinct and less 
explored regime in which the system qubit resonantly couples strongly to only a few TLS 
defects, as opposed to a dense spectrum. Similarly, it may be the case that, given the much 
smaller volumes (~ 1CT 18 m 3 ) of the micronscale mechanical resonators currently under in- 
vestigation, the relevant system-environment model is an oscillator interacting with only a 
few TLS's. In the following sections we shall analyze just such a system. 

In Sec. [II] we give some simple estimates based on existing bulk system TLS theory in 
order to motivate our mechanical resonator mode-few TLS model as well as to anticipate 



some of the consequences of TLS-dominated mechanical damping/decoherence. In Sec. Ill 



we derive the model closed system resonator-TLS Hamiltonian, and in Sec. |IV| we present the 
open system master equation with further details of the derivation given in the Appendix. 



Section M focuses on the effect of a few TLS's on the resonator damping, while Sec. VI 



describes the consequences for the decoherence of mechanical resonator superposition states. 



Section VII provides a few concluding remarks. 

Although the present paper focuses on the role of TLS's for micronscale mechanical oscil- 
lator damping/decoherence, we do not completely neglect other mechanisms. In the spirit 
of keeping our model as simple as possible, we lump together all other relevant damping 
and decoherence mechanisms, such as clamping loss j 28 l 29 l 3 ° l 31 l as an additional oscillator bath 
to which the mechanical system mode couples. This will allow us to gauge somewhat the 
extent to which other baths 'interfere' with the TLS bath in their damping and decoher- 
ence effects on the oscillator system. For example, the system oscillator's net damping rate 
need not be the sum of the damping rates due to the individual baths. One point that 
should be emphasized in this context is the highly nonlinear, quantum nature of the coupled 
oscillator-TLS (equivalently spin-1/2) dynamics. Exact analytical or even simpler approxi- 
mate equations are hard to come by and so we will resort to solving for the full dynamics 
using numerical methods. We will be limited computationally to considering only a few 



3 



TLS's-three to be precise. In future work we plan to find ways to analyze the effects on 
mechanical damping/decoherence of larger numbers of TLS's. 

Another source of nanomechanical damping and decoherence that we do not explicitly 
take into account is the measurement process itself J22E3EI] The resonator damping described 
in Sec. [V] can be probed using, for example, continuous in time position detection with a 
single electron transistor We shall assume that the resulting back reaction on the resonator 
due to the position detector can be simply modeled by the same additional oscillator bath 
at some finite temperature.^^ The resonator superposition state decoherence described 



in Sec. VI can be probed using, for example, the microwave cavity-superconducting qubit 



scheme outlined in Refs. 

In the present paper we neglect mechanical strain (i.e., phonon) mediated coupling be- 
tween TLS's, assuming the latter to couple directly only to the oscillator system mode 
and with the TLS's damping treated phenomenologically, characterized by a decay time T\. 
Given our current, almost complete lack of theoretical understanding of the role of TLS's 
for damping and decoherence in nano-to-mesoscale mechanical resonators, we feel that it 
is worthwhile to start with this simpler, noninteracting TLS model. The low temperature 
acoustic pulse probe investigations of bulk amorphous solids^^ and the mechanical qual- 
ity factor and resonant frequency measurements of much larger resonator d 39 * 40 ! point to the 
likely importance of interactions between TLS'sP^ 42 * 43 ! For example, phonon echo experi- 
ments yield TLS dephasing times that are much shorter than their T\ lifetimes, thought 
to be due to TLS spectral diffusion arising from the phonon mediated interaction between 
non-resonant TLS's. 41 With the significantly reduced volumes of nano-to-mesoscale mechan- 
ical resonators, these strain interactions may in fact be considerably enhanced. We plan to 
analyze the effects of such TLS interactions on nanomechanical damping/decoherence in a 
future work. 



II. SOME ESTIMATES 

In this section we adopt various results from earlier analyses of bulk, amorphous systems 
to try to gain some initial idea of expected consequences for damping/decoherence of nano- 
to-mesoscale mechanical resonators due to TLS's. We begin by estimating the number 
magnitude of TLS's that are near resonance with the mechanical mode frequency of interest, 
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which we shall in this paper assume to be the lowest, flexural mode. For a range of bulk, 
amorphous solids, experiments are consistent with a TLS distribution of the forrrJ-^ 

dN = V^-dA dA b , (1) 

where V is the mechanical resonator volume, Ao and are the asymmetry and tunnel 
splitting energies of the TLS's potential double well (see Fig. JTJ, and P is the approximately 
constant spectral density that can be expressed as 

P = C^. (2) 

Here, p is the mass density, v is the deformation potential (approximated as isotropic) and v 
is the speed of sound (approximated as isotropic and polarization independent). The dimen- 
sionless constant C ~ 10~ 4 — 1CT 3 is approximately universal J 44 * 45 ! Of course, the nanome- 
chanical resonator may not be fashioned out of one of the amorphous materials surveyed 
in Ref. [15], but instead out of a crystalline material. In such a case we view Eq. Q and 
the estimates we shall now derive for the total TLS number as providing an upper bound. 
Integrating Eq. Q to obtain the total number of TLS's in energy width 5E about the TLS 
eigenenergy E = a/Aq + A 2 , we obtain: 

■i + y i-( Wg) . i 

1 - VI - (SW-E) 2 ] 

where E mm is the low energy cutoff in the TLS distribution.^ Taking the typical ball- 
park values pv 2 ~ 10 11 kgm _1 s -2 and u ~ 1 eV gives for the spectral density ([2]), 
P ~ 10 44 -10 45 J-^" 3 . Substituting this into Eq. Q and taking E = fko [i.e., m resonance 
with the oscillator mode u), 5E = Huj/Q, where Q is the oscillator mode quality factor, 
and E m [ n /kB = 1 mK, we obtain for our estimated total TLS number close to resonance 
expressed in natural units: 

8N~ 0.1-1 (^ 3 )(//GHz) 

(Q/10 4 ) 1 ' 

Thus, according to this estimate, micron-scale, radio frequency mechanical resonators are 
on the borderline between being unlikely to have a single TLS close to resonance and being 
very likely with a moderate increase in size. 

We next estimate the nanomechanical fundamental mode displacement amplitudes that 
saturate the TLS's. Observation of quantum effects in the mechanical dynamics will require 



SN = -VP In 
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cooling the resonator such that ksT < huj. In this regime resonant absorption damping 
is expected to dominate. Theory of resonant TLS mechanical damping and acoustic pulse 
decay in large amorphous resonators and bulk solids gives for the energy damping rate^ 

_, ixCuj ( hw \ 

r damp « 7= ; , ■ t&nh — — , (5) 



where C is the same dimensionless constant as introduced above, e is the elastic strain 
amplitude (approximated as isotropic), v is the deformation potential, T x is the TLS damping 
time, and T 2 is the TLS transverse relaxation time, related to the dephasing time as 
T 2 X = (2T1)- 1 + T^ 1 . From this expression we see that the saturation threshold strain is 
given by 

(6) 



Vy/TiT 2 

For the example of a doubly-clamped beam mechanical resonator of length I, thickness d, 



(Y/A) ~ 10- 6 [L //i m } (7) 

'■'/(N) (d/fim) y (T\T 2 / )isec 2 ) 



and midpoint 1/2 transverse displacement amplitude Y (see Sec. Ill), the volume averaged, 
rms strain is e ~ dY/l 2 . Substituting this into Eq. ^ and expressing in terms of natural 
units, we obtain 

(/7/im 2 ^ 

Earlier phonon echo experiments in fused silica glass found TLS damping and transverse 
relaxation times T\ ~ 100 /isec and T 2 ~ 10 /isec, respectively, at T w 20 mXP Thus, 
the transverse relaxation is dominated by dephasing: T 2 w T^. Given that quantum zero- 
point displacement uncertainties of micron-scale mechanical resonators are ~ io- 4 -io~ 3 AP 
which, as can be seen from Eq. (|7]), exceed the saturation threshold, we thus expect that 
experiments which measure damping and decoherence of such resonators will operate well 
into the saturation regime. As can be seen from Eq. ([5]), one important consequence is that 
resonant TLS dominated damping is expected to be amplitude dependent, corresponding 
to a nonlinear damping force that depends on both the position and velocity coordinates of 
the mechanical oscillator. 

Our final approximation concerns the decoherence of mechanical superposition states. 
For weak mechanical damping (i.e., Q ^> 1) and provided temperatures are not too low (i.e., 
ksT S> Hlu/Q), the effective dynamics of an oscillator bilinearly interacting with a bath 
of free oscillators satisfies a quantum fluctuation-dissipation relation between the system 
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oscillator's damping and decoherence rates: 



1 / AY\ 2 , / huj 



Tdecohere — ^damp ^ ~y~ J COt ^ (^A^T/ ' ^ 

where AY is the uncertainty in the oscillator's position and Y zp = y/h/ (2mw) is the oscil- 
lator's quantum zeropoint position uncertainty. Assuming that this fluctuation-dissipation 
relation applies also to the bulk oscillator with the TLS bath, we can substitute in expres- 
sion ([5]) for the damping rate to obtain the decoherence rate. The first thing to notice is 
the cancellation of the hyperbolic temperature functions. Below the saturation threshold we 
would therefore conclude that the decoherence rate is temperature independent while the 
damping rate is temperature dependent for an oscillator coupled to a TLS bath. In this 
respect, the oscillator system-TLS bath is 'dual' to the oscillator system-oscillator bath in 
the sense that the decoherence rate of the latter has the inverse temperature dependence of 
the damping rate of the former, while the damping rate of the latter and decoherence rate of 
the former are both temperature-independent.^ However, this duality is to a certain extent 
academic since the above saturation estimates suggest that nano-to-mesoscale mechanical 
resonators will be well within the saturation regime for quantum superpositions of distinct 
position states that must necessarily be larger than the zeropoint position uncertainty. As 
a consequence, the temperature dependences of the TLS T\ and T 2 dependent terms must 
also be taken into account. Since a good understanding of the relaxation mechanisms of 
TLS's is lacking in nano-to-mesoscale mechanical resonators,^ we will not attempt to make 
predictions for the temperature dependences of various observable quantities in the present 
paper. 

Leaving aside temperature dependencies, another notable consequence of applying the 
above quantum fluctuation-dissipation relation to the oscillator system-TLS bath is the 
weaker (i.e., linear) dependence of the decoherence rate on oscillator position uncertainty. 
Thus, at low temperatures we might expect the decoherence rate to increase more gradually 
as the position separation in the quantum superposition state is increased, as compared with 
the quadratic separation dependence for the oscillator bath. Of course, other decoherence 
mechanisms, e.g, due to clamping loss, will then be expected to eventually dominate if they 
have the stronger quadratic dependence. 

Summarizing the findings of this section, we expect that a relevant model for a nano- 
to-micronscale mechanical resonator interacting with TLS's is a single oscillator coupled 



to a bath of a few TLS's. Furthermore, we expect that mechanical damping will show an 
amplitude dependence, while the decoherence of mechanical superposition states will depend 
weakly on their position separation. The following sections will bear out these expectations. 

III. RESONATOR-TLS HAMILTONIAN 

In this section we will derive the Hamiltonian describing the dynamics of the lowest, 
fundamental flexural mode of a doubly-clamped beam mechanical resonator interacting with 
TLS's that are located randomly throughout the beam volume. Related analyses are given 
in Refs. |19|50] . We shall assume a long, thin elastically isotropic beam with length I, width 
w, and thickness d satisfying / ^> w > d, mass density p, and bulk modulus K (Fig. [T]). The 




FIG. 1: Flexing doubly-clamped beam (exaggerated for clarity). Inset: TLS double well potential. 



equation of motion for small transverse displacements y(x,t), < x < I, of the beam i 

where b = d/^/l2 is the bending moment and we assume zero applied longitudinal strain. 
The total energy of the beam is 

E = -pwd [ dx{yf + -Kwdb 2 [ dx(y"f. (10) 
2 Jo 2 y 



Solving Eq. ^ with clamped boundary conditions y(0) 
obtain for the lowest frequency (fundamental) eigenmode: 



y(x,t) = Y(t) 



4>{x) 



y(l) = y'(0) = y'{l) = 0, we 



(11) 



where the normalised eigenfunction is 



x) = A | sm ^ — - — j — smh 



/nax 



+ (3 



/irax\ /7rax\ 
(■( >s I — — ) — cosh 



I 



(12) 



with a ~ 1.51 obtained from the clamped boundary condition expression cos(7ra) cosh(7ro;) = 
1 and f3 = [cos(7ra) — cosh(7ro;)]/[sin(7ra) + sinh(7ra)] w —1.02. The constant A pa 0.983 is 
fixed by requiring that 4>(x) be normalized as follows: 



n 1 / dx(p 2 (x) = 1. 
Jo 



(13) 



The time-dependent part of the solution (11) is Y(t) = Re(Yoe lujt ) , where the fundamental 
mode frequency is 



{na\ 2 IK 



6.46 



d K 

^V7' 



(14) 



The solution (11) is expressed such that Y(t) gives the transverse displacement of the beam 



at its midpoint x = 1/2. Substituting Eq. (11) into the total energy (10) and employing the 



normalization condition (13), we obtain 

1 



E = -m(j)- 2 (l/2)Y 2 + -m<J)- 2 (l/2)uj 2 Y 2 , 

_ _ 



(15) 



where m = pwdl is the mass of the resonator. Thus, the fundamental mode dynamics is 
that of a harmonic oscillator with effective mass 



m m 



2 (//2) 2.52 : 



(16) 



where, from Eq. (12), we have used 4>(l/2) ~ 1.59. 



Quantizing the fundamental mode, we introduce mode raising and lowering operators a + , 
a; [a, d + ] = 1, with 

Y = Y zp (a + a + ), (17) 

where Y zp = y/h/ (2m c guj) is the quantum zeropoint displacement uncertainty. The free 
beam, fundamental mode Hamiltonian is then simply H m = hw(d + d + 1/2). 
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Moving now to the TLS Hamiltonian, we have: 



N 

E 



\^of + Ufa® 



(18) 



where j = 1, 2, iV labels the TLS, Aq is the asymmetry of the j'th TLS's potential well 
and A® is its tunnel splitting that depends on the well barrier height and width (see Fig. [lj). 
Mechanical resonator motion couples to a TLS largely through the strain dependence of the 
asymmetry energy: 

(19) 



,00 

z kl ) 



fc,Z=l 

where z/^/ is the deformation potential and eu is the elastic strain tensor. For small amplitude, 
transverse y(x, t) flexural displacements of long, thin beams, the nonvanishing strain tensor 
components for a defect located at x and a distance — d/2 < £ < d/2 normal to the neutral 
(i.e., strain-free) surface (see Fig. IT]) are e xx ~ —(d 2 y/dx 2 and e yy = e zz ~ a(d 2 y/dx 2 , where 



a is Poisson's ratio. From Eqs. (j X 1 [) and (|17|), the transverse displacement field operator is 

(20) 



4>(l/2) 



Subsituting (20) into (19), we obtain for the mechanical resonator-TLS defect interaction 
Hamiltonian 

N 

- (21) 



Ant = ^A^(a + a+)aW, 

where the resonator-TLS strain coupling strength for defect j located at (xj, (j) in the beam 
takes the form 



\V) = -Y zp y{\ - 2a)Q 



(fix,) 



(22) 



~"'- J <j>{l/2) 

and we have assumed for simplicity an isotropic deformation potential coupling v. Finally, 
writing out the full resonator-TLS system Hamiltonian, we have 



N n 1 

£ 5 = hw(d + d + 1/2) + ^ -A?Vjp> + -A? V« + A^(g + G + )a^ 



i=i 



(23) 



The strength of the coupling A 1 -^ depends on the location of the TLS defect. In particular, 
the coupling is strongest for a defect on the surface at the beam ends, i.e., Q = ±d/2 and 
Xj = 0, 1. In order to gain a sense of the expected magnitudes of the coupling, it is convenient 
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to express the various beam material constants and dimensions in natural units. We obtain 
for the dimensionless coupling strength: 

^ . 10 _3 fb\ fr{*i)\ ( p V'V K Y m (< l - 2 «) 



tku \d J V 0"(O) J V.10 3 kgm"V V10 n Nm~V V eV 

J±Y*(^(±Y\ (24) 



Jim. J \/Um/ \^ m , 
while the fundamental flexural mode frequency / = u/(2ir) expressed in natural units is: 

(gL) W L3 (l0 3 kgm- 3 ) (lO^Nm- 2 ) (/xm) (^m) ' (25) 
IV. OPEN SYSTEM MASTER EQUATION 

In this section we derive the master equation for the coupled resonator-TLS system, taking 
into account the environment of the system. In an actual beam mechanical resonator, the 
fundamental flexural mode will couple not only to the TLS defects, but also to the other, 
higher frequency resonator modes via anharmonic interaction terms. The fundamental mode 
will also couple to bulk, substrate modes at the beam supports. Furthermore, the TLS's will 
couple to the higher frequency resonator modes through the strain dependence of the TLS's 
asymmetry energies. The latter will not only cause damping of the TLS'sj^ but will also 
induce interactions between the TLS'sP^* 43 * 52 ! However, our goal in the present investigation 
is not to accurately model the respective environments of the fundamental mode and TLS's, 
but rather as a first step to consider the simplest possible idealized model environments in 
order to gain an idea of the quantum dissipation and decoherence dynamics of the mechanical 
resonator interacting with damped TLS's. 

As idealised model environments, we consider baths of non-interacting harmonic oscilla- 



tors. Hamiltonian (23) is then augmented by the environment Hamiltonian and coupling 
term: 

JV 

H cnv + Hs-cnv = ^nbn + + Y ^ K n(K + &n) 

n j=l n n 

N 

+E^' ) E^ i) ( c - )+ + c « ) )' ( 26 ) 

j=l n 

where recall Y = Y zp (a + a + ) and for notational convenience we have dropped the hats 
on the operators and have also neglected the bath zeropoint contributions. In our model 
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each TLS is assumed to couple with strength Kft to independent, noninteracting oscillator 
baths, characterized by environment mode operators Cn\j = 1,2,...,N. The environment 
mode operators b n couple directly to the resonator with strength K n , collectively modeling 
all energy loss mechanisms other than those involving the TLS's, such as clamping loss and 



anharmonic processes. Combining Hamiltonians (23) and (26), we have the total system- 
environment Hamiltonian H = Hs + H env + Hs-env 

In Appendix [A] we apply the self-consistent Born approximation together with a Markov 
approximation to obtain the following master equation describing the dissipative dynamics 
of the coupled resonator-TLS system: 

Ps(t) = -ll H s^s(t)]-g[ydPy,Ps(t)}]-^coth^^^ [Y,[Y, Ps (t)]] 
A - / /?(i)\ 2 



N 1 



-£^)(5) tanh (i^)^'K«.^)}], P7) 

where psif) is the resonator-TLS system density matrix, Py is the resonator momentum, 
{■, ■} denotes the anticommutator and = \J (A^) 2 + (A^) 2 is the jth TLS energy level 
separation. The parameter 7 gives the energy damping rate of the resonator in the absence 
of the TLS, while gives the jth TLS relaxation time from its excited energy eigenstate 
in the absence of the resonator. 



One potential advantage of the damping time/rate parametrization used in Eq. (27) is 
that it does not in fact depend explicitly on the microscopic nature of the environment and 
how it couples to the system. As long as the Markov approximation can be made and the 
fluctuation-dissipation theorem holds, then the damping and diffusion terms are uniquely 



related so that Eq. (27) can be assumed to apply for other system-environment interactions 
as well. This then allows environment model-independent predictions for the open system 
dynamics provided that they are expressed in terms of the damping parameters 7 and , 
as opposed to predictions concerning the explicit temperature dependence,^ which depend 
on the nature of the TLS environment. 
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V. DAMPING 



In this and the following section, we present the results of numerically solving the master 
equation (27) using the Quantum Optics ToolboxP^Most of the results are for the mechanical 
resonator mode subsystem only with the TLS sector of state space traced over; we assume 
that it is the resonator mode which is directly probed in experiment. 

We begin this section with a focus on the damping of the mechanical resonator coupled to 
a single TLS. Both the resonator and the TLS are coupled to independent ohmic oscillator 
baths. We assume that the system is initially in a product state, ps(0) = p res (0) ® Ptls(O), 
where the TLS is initially in a thermal state Ptls(O) = exp[— if TLS /fc B T]/Z with .Htls 



defined in Eq. (18) and Z = Tr (exp [ — ./^tls /k bT\ ) . Similarly, the initial resonator state is 
a thermal state with H m = fkd(aJa + 1/2) that has been displaced using the operator d = 
exp[xo(a^ — a)/2], where xq is the initial displacement of the thermal state from equilibrium 



in units of the quantum zero-point displacement uncertainty Y zp defined in Eq. (17); from 
now on we use x to denote the mechanical resonator fundamental mode displacement in 
units of the zero-point uncertainty. 



We use Eqs. (24) and (25) to determine the resonator-TLS coupling constant A and the 
fundamental flexural mode frequency /. For I = 10 pm, w = 1 pm, d = 0.1 pm, we find 
X/tvjj « 10~ 2 and / « 100 MHz, or u « 2ir x 100 MHz. Where convenient, we shall use 
dimensionless time units, t — > ut, with T\ and 7 expressed as uTi and j/uj, respectively, 
and A, A i; and temperature T expressed as X/hui, Ai/hw, and ksT jhuj, respectively. 

We begin with a symmetric TLS (Ao = 0) that is on resonance with the mechanical 
oscillator (A& = 1). In Fig. [2j we plot the envelope of the resonator's ensamble averaged 
position versus time, corresponding to the so-called interaction picture or rotating frame, 
where the resonator and TLS's rapid free evolution are factored out. We give the resonator 
an initial displacement Xq = 3 and assume a range of experimentally realistic values for the 
TLS relaxation time that are relevant at mK temperatures, 1 psec < T\ < 100 psecp^ and 
also assume a non-TLS resonator energy damping rate 7 = 10~ 5 . Comparing the envelope 
of the resonator's motion at T = 10 mK (= 2.02 in dimensionless units) for T\ = 100 psec 
(= 6.5 x 10 4 in dimensionless units) (curve B) to the envelope of the resonator in the absence 
of the TLS (curve A), we see that the TLS and its bath cause significant amplitude damping 
of the resonator, even for large values of T\. Fig. [2] also shows an increase in the damping of 
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the resonator for decreasing values of Ti; for shorter TLS damping times, energy exchanged 
between the resonator and the TLS is more quickly dissipated through the TLS bath. 




FIG. 2: Envelope of the motion of the resonator coupled to an on resonance TLS for T\ ~ 6.5 x 10 4 
(B), T\ ~ 1.6 x 10 4 (C), and T\ ~ 6.5 x 10 2 (D). The solid black curve (A) shows the damping of the 
resonator with no TLS coupling. For all curves A = 10~ 2 , Ao = 0, A& = 1, 7 = 10 -5 , T = 2.02. 

Note also from Fig. [2] that the resonator's average position damps out and then revives 
on time scales much less than the TLS relaxation time, Ti, appearing to indicate a complete 
transfer of energy from the mechanical resonator to the TLS. However, the on-resonance 
TLS can only absorb a maximum of one quantum (phonon) of vibrational energy, while the 
energy stored in the resonator for the considered initial displacement Xq = 3 corresponds to 
more than one phonon on average. This apparent paradox is resolved by noting that the 
quantum resonator's average energy depends on the average of the position squared, not 
the square of the average; not all of the energy is transferred to the TLS when the average 
position completely damps out. This can be seen in Fig. [3j where the individual resonator 
and TLS energies, as well as the total TLS+resonator energy, are plotted over time. 

Energy is able to be transfered between the oscillator and the TLS because of the long 
TLS T\ times. If T\ is shorter then the energy transferred to the TLS dissipates to its bath 
more rapidly, with the result that less energy returns to the oscillator. The effect of a shorter 
T\ on the oscillator and oscillator+TLS energy is shown in Fig. |3j Fig. [4] demonstrates the 
dependence of the resonator dynamics on the oscillator-TLS coupling constant A. For weaker 

14 



couplings the damping out and revival of the oscillator's motion occurs at later times. 




FIG. 3: Dependence of resonator and TLS average energy vs. time for the same parameters as 
those used in Fig. [jjj with T x « 6.5 x 10 4 (black) and T x xs 6.5 x 10 2 ( green). (A, D) Resonator 
energy; (B, E) TLS energy; (C, F) Resonator + TLS energy. 



Because a TLS absorbs only a finite amount of energy E = a/Aq + A 2 , it can become 
saturated and hence affect the resonator damping rate. In Fig. [5] we plot the resonator 
amplitude decay rate 7a versus time for a range of initial displacements xq, where ja is 
defined as follows: 

7A = -^yln((xM))). (28) 

In order to display trends in damping more clearly we have increased the values for 7, A and 
T^ 1 , shortening the time over which the damping takes place. We see in Fig. [5] that as Xq 
increases, the initial decay rate of the resonator decreases. For a large initial displacement, 
xq = 10 (D), the initial decay rate approaches that of the resonator in the absence of the 
TLS (E), suggesting near total saturation of the TLS. For later times, however, when the 
resonator amplitude has decayed to a value near the zero point displacement, the decay 
rates for all four initial displacements approach a common, constant value. For xq = 5 this 
value is reached at out ~ 250, while for xq = 10 this value is reached at ut w 375. We 
know that the oscillator bath's contribution to damping is constant at all times (Fig. [5jE]), 
whereas the resonator coupled to the TLS shows amplitude-dependent damping initially, and 
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FIG. 4: Dependence of the resonator ensemble-averaged amplitude vs time for different resonator- 
TLS coupling constants A. The parameters are the same as those used in Fig. [2] with T\ ~ 6.5 x 10 4 , 
A = 0.125 (A), A = 0.1 (B), A=0.075 (C), and A = 0.05 (D). 
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FIG. 5: Left: Amplitude decay rate of resonator coupled to TLS with T\ = 10 for xq = 3 (A), xq = 5 
(B), xo = 8 (C), xq = 10 (D) and without TLS (E) as a function of time. Right: Dependence 
of crossover amplitude xco on T\ for xq = 8. For all curves A = 0.1, Ao = 0, A5 = 1, 7 = 
0.01, andT = 1. 

amplitude-independent damping at later times (Fig. [5|A-D]). These two distinct behaviors 
indicate that the TLS is indeed saturated at higher resonator amplitudes, while for lower 
amplitudes the TLS is unsaturated and its contribution to damping is therefore uniform. 
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On the right-hand side of Fig. [5] we investigate the T\ dependence of the amplitude at which 
the crossover from amplitude-dependent to amplitude-independent damping occurs. The 
crossover amplitude xqq is plotted as a function of 7\, showing a nearly linear relationship. 

Fig. [6] shows the damping of a resonator with initial displacement xo = 5 coupled to 
both an oscillator bath and a TLS (A), and coupled to the TLS (B) and oscillator bath 
(C) individually. Curve D is the sum of curves B and C. The substantial difference between 
curves A and D demonstrates that one cannot simply add the individual damping rates to 
obtain the net TLS+oscillator bath damping rate when both these sources are present. 

0.04 
0.035 
0.03 
0.025 
f< 0.02 
0.015 
0.01 
0.005 


50 100 150 200 250 300 350 400 

FIG. 6: Decay rate of resonator for xq = 5 with TLS and oscillator bath (A), with TLS only 
(B), with oscillator bath only (C), and sum of curves B and C (D). For all curves A = 0.1, Ao = 
0, A b = 1, 7 = 0.01, T = 1, and T x = 10. 

So far we have considered just the special case of a symmetric (A = 0) and resonant 
(A;, = 1) TLS. In an actual mechanical resonator there will be a distribution of TLS's, 
not necessarily symmetric or on resonance [see Eq. (JTJ]. In Fig. [7J we show the resonator's 
amplitude damping rate 7a for a range of TLS Ao and A& values. In the absence of a TLS, 
7a is equal to half the energy damping rate, 7, due to the oscillator bath. Fig. [7] shows 
that 7a is greatest for A = and A&=1, and sharply decreases to 7/2 as A^ is moved off 
resonance. The resonator's damping rate decreases more gradually as the TLS is made more 
asymmetric (A 7^ 0). The damping rates in this figure were extracted from an oscillator 
with small initial displacements in the range of constant values shown in Fig. |5j and therefore 
showed no time or amplitude dependence. 

We now consider the damping of the resonator coupled to three TLS's, which have energies 
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FIG. 7: Resonator damping rate for various values of Ao and The parameter values A, 7, T, 
and T\ are the same as those used in Fig. [5j 

close to the resonator energy. The TLS's A3 values are indicated by the diamond, circle, and 
triangle symbols in Fig. [8} The TLS energies were selected randomly from within the range 

0. 75 < E < 1.25 using the distribution ([!]). As is the case for a single TLS, we expect the 
damping of the resonator when coupled to multiple TLS's to show amplitude dependence 
for early times due to TLS saturation. Fig. [9] shows the decay in the presence of the three 
TLS's and the oscillator bath for three different initial displacements x$ (curves A, B, C). 
As expected, the resonator decays more quickly for smaller initial displacements, while for 
larger initial displacements the decay rate approaches that due to the oscillator bath only, 
indicating saturation. Unexpectedly, at later times there is a crossover to a much slower, 
constant decay rate that is even less than that of the resonator in the absence of the TLS's, 

1. e., due to the oscillator bath only. The unnormalized amplitudes of the envelope at which 
the crossover occurs are the same for each curve, independent of the initial Xq displacement. 
We have also found similar crossover behavior in amplitude decay rates assuming other 
randomly selected distributions of three TLS's with energies close to the resonator energy. 
The bumps in curves B and C are a result of numerical approximation and are not a product 
of the system's behavior. 
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FIG. 8: Considered distribution of TLS's. The energies of the TLS's are Ao = 1.2002 and Aj, = 
0.1475 (diamond), A = 0.6108 and A b = .9145 (circle), A = 0.0070 and A b = 1.0354 (triangle). 
The dashed line denotes the TLS energy on resonance with the mechanical oscillator. 




FIG. 9: Damping of the oscillator coupled to three TLS's with initial displacements xq = 3 (A), 
x = 5 (B), and x = 8 (C). For all curves A = 0.1, 7 = 0.01, T = 1, and T x = 10. (D) shows the 
damping of the oscillator in the absence of the TLS's. 
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VI. DECOHERENCE 



In this section we investigate the contribution of the TLS's to the decoherence of a 
mechanical resonator Schrddinger cat state. The resonator's initial state consists of a 
superposition of two coherent states: \ip) = N' 1 (| + xq) + \ — Xq)), where | ± Xq) = 
exp(±xo(a^ — a)/2)\0) with initial displacement ±xo and iV is a normalization factor. The 
parameter choices are 7 = 0.01, A = 0.1 and T\ = 10, the same as those used for showing 
the trends in damping in Sec. |V| In order to study the evolution of the decoherence of the 



cat state, we evaluate the Wigner function of the resonator density matrix p. Fig. 10 shows 
the Wigner function in the absence of TLS's, i.e., with only the ohmic oscillator bath, at 
five equal time step increments. Within a single period of the oscillator's motion, the in- 
terference fringes between the two cat states decay substantially for the assumed parameter 



choices. Fig. 11 shows the evolution of the Wigner function for the same initial resonator 
state but with a single, on-resonance TLS coupled to the resonator in addition to the ohmic 
bath. The interference fringes decay more rapidly due to the presence of the TLS. Fig. 12 
shows the Wigner function for the resonator coupled to a single TLS, but in the absence of 
the ohmic oscillator bath. The decoherence time in this case is longer than the resonator's 
oscillation period. 
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FIG. 10: Wigner function of the resonator density matrix in the absence of the TLS, where 
x = 5, 7 = 0.01, T = 1, and Ty = 10. 



In order to obtain a more quantitative understanding of the TLS-induced decoherence, 
we take the average amplitude of the interference fringes in a small, disk-like region centered 
between the two peaks, and plot the negative log of this amplitude, which we denote yl int , 
as a function of time. This quantity is essentially the same as that used in Ref. [HI], where 
they show that the decoherence rate of a harmonic oscillator cat state coupled to an ohmic 



oscillator bath is proportional to the slope of A int . Fig. 13 shows A int as a function of 



20 





'. / .-J 




iut-.T'4 








Mil— j.t;4 




-_ Y 

'. .'I-.T 


6 




t 




I 




1 




E 




2 


# 


■1 

2 


w 


4 

z 




4 

2 


■ 


2 




« 3 


# 






x {I 




* £ 




x 




-2 


• 


-2 
-4 


• 


-2 
-4 




-2 
-4 




-2 

-4 




-'6 


-k 




-6 




-fi 




-6 


1 # " 




-S -4 -2 14 6 
P 




-6 -4 -2 2 4 6 
P 




-6 -4 -2 o a 4 6 
P 




-4-4-2 £ i 6 
f 




-fi -4 -2 2 4 6 
P 


FIG. 


11: Wigner function for the resonator coupled to the TLS, where xq = 
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FIG. 12: Wigner function for the resonator coupled to the TLS, but without the resonator's 
independent ohmic oscillator bath, where xq = 5, A = 0.1, Aq = 0, A& = 1, T = 1, and T\ = 10. 



Note: time steps are larger than for Figs. 10 and 11 




FIG. 13: j4i nt vs. ut for the resonator in the absence of TLS's. Initial displacements are xq = 5 (A), 
x = 5.5 (B), x = 6 (C), xq = 6.5 (D), x = 7 (E), with A = 0, A b = 1, 7 = 0.01, T = 
1, and T% = 10. 
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FIG. 14: A[ nt vs. cot for the resonator coupled to a single, on-resonance TLS, but not to the 
resonator's independent ohmic oscillator bath. Initial displacements are xo = 5 (A), xo = 5.5 (B), 
x = 6 (C), x = 6.5 (D), x = 7 (E), with A = 0.1, A = 0, A b = 1, T = 1, and T x = 10. 

time for the resonator in the absence of the TLS. As predicted by theory, we can extract 
a decoherence rate from the constant increase of A int {uot). We find that the decoherence 
rate, or the slope of A int (ujt), goes as the square of the initial displacement, as expectedP 



Fig. [14] shows A[ nt as a function of time for the resonator interacting with an on-resonance 
TLS, but without the resonator's independent ohmic oscillator bath. The decoherence due 
to the TLS does not show the same dependence on the initial displacement as the oscillator 
bath-induced decoherence; in fact, there is no apparent systematic dependence on initial 
displacement. In particular, increasing the initial displacement does not result in an overall 
increase in the decoherence rate. 



In the Wigner function plot Fig. 12 for the resonator coupled to the damped TLS only, 



the interference fringes in the center region partially decay away and then return after a 



full mechanical period. These interference oscillations, which can also be seen in Fig. 13 
are a consequence of the on-resonance TLS behaving like a position measuring device. As 
the resonator interacts with the TLS, the two distinct position states in the superposition 
cause the TLS state to evolve in different ways, resulting in an entangled resonator-TLS 
state; the resonator partially decoheres. However, because the coupling between the TLS 
and its environment as parametrized by the T\ time is relatively weak, subsequent evolution 



22 



almost completely undoes the entanglement; the resonator partially recoheres, restoring the 
interference between the two position states. 



In Fig. 15 we plot the decoherence of the resonator coupled to three TLS's for three initial 
displacements (A, C, and E), as well as the decoherence of the resonator coupled only to 
the oscillator bath for comparison. The decoherence due to the presence of the three TLS's 
is greater than that due to a single TLS, while the decoherence due to the oscillator bath 
shows a systematic dependence on displacement that the TLS-induced decoherence does not 
exhibit. Decoherences and partial recoherences on the shorter, mechanical period timescales 
can be clearly seen, due to the mechanical resonator-TLS's entangling- disentangling dy- 
namics as discussed above for the single TLS case. Similar behavior was found for other 
randomly selected sets of three TLS's with energies near the resonator energy. Finally, in 



Fig. 16 we plot A mt vs. ujt for the resonator coupled only to the ohmic bath (A), coupled 



only to the three TLS's (B), and coupled to both the ohmic bath and the three TLS's (C). 
As expected, we find that decoherence occurs more rapidly when the oscillator is coupled to 
both the ohmic bath and to the three TLS's. 




FIG. 15: Displacement dependence of A^ nt vs. oot for the resonator coupled to the independent 
ohmic oscillator bath (A, C, and E) and to the three TLS's with energies indicated in Fig. [8] without 
the oscillator's ohmic bath (B,D, and F). The initial displacements are xq = 5 (A, B), xq = 6 (C, 
D), and xq = 7 (E, F). The parameter values are A = 0.1, Aq = 0, A& = 1, T = 1, and T\ = 10. 
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FIG. 16: A[ Qt vs. tut for the resonator coupled to the independent ohmic oscillator bath only (A), 
to three TLS's only (B), and to the ohmic oscillator bath and the three TLS's (C). The TLS Aj 



values are the same as those in Fig. 15 The initial displacement is xq = 5, and the parameter 



values are A = 0.1, 7 = 0.01, T = 1, and T\ = 10. 
VII. CONCLUSION 

In the present paper we have investigated the damping and decoherence dynamics of the 
flexural mode of a nanomechanical beam resonator interacting with a few damped TLS's. 
It was found that the resulting damping rate is amplitude dependent, while the decoherence 
rate for superpositions of position states depends only weakly on their separation. This is to 
be contrasted with the damping and decoherence trends of the more commonly considered 
resonator interacting bilinearly with an ohmic bath of free oscillators. In the latter case, 
the resulting damping is amplitude- independent, while the decoherence rate scales with the 
square of the position separation in the initial superposition state. 

In our model, strain-mediated interactions between the TLS's were neglected. It will 
be interesting to take such interactions into account; the reduced volume of a suspended, 
nanoscale resonator may result in a significant enhancement of the strain-mediated TLS-TLS 
interaction as compared with the TLS-TLS interaction in the bulk. It is also important to 
try to increase the number of TLS's in our numerical experiments and compare the obtained 
damping and decoherence trends with approximate analytic results derived for the model 
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system of a resonator interacting weakly with a dense spectrum of TLS's. This will then 
enable a test of the analytical approximations, as well as help shed light on the role of TLS's 
in the larger, optomechanical resonator experiments. 
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APPENDIX A: DERIVATION OF THE OSCILLATOR-TLS DISSIPATIVE MAS- 
TER EQUATION 

In terms of the system-environment density matrix p and Liouvillian superoperator L 
defined by Lp = [H, p] where H is the total system-environment Hamiltonian, the time- 
dependent Schrodinger equation is 

p = ~l p . (Al) 

We wish to derive a master equation for the system density matrix comprising the oscillator 
and TLS only: p$ = Tr env p, where the trace is performed over the oscillator-TLS environ- 
ment. Following the Nakajima-Zwanzig projection operator methocP^ along the lines of 
Ref. [57], we introduce projectors 

P(-) = Pcnv ® Tr cnv (-) (A2) 

and 

Q(-) = (l-Pcnv®Tr cnv )(-), (A3) 

where p env = e~^ Hcnv /Z , with (3 = 1/ksT, for environment temperature T. Suppose that at 
the initial time t = 0, the system and environment are in a product state p(0) = ps(0) ®p env . 
Then Pp(0) = p(0) and Qp(0) = 0. Now partition the density matrix p and Schrodinger 
equation using the projectors P and Q: 

p x = - l -PL( Pl + p 2 ) (A4) 
P2 = - l -(l-P)L( Pl +p 2 ), (A5) 

25 



where p\ = Pp and P2 = Qp. Solving formally for p 2 , we have 



P2(t) = / dt'exp 



[l-P)L Pl (t-t') 



(A6) 



where we have used the fact that pzfi) = Qp(0) = 0. Substituting Eq. (A6) into Eq. (A4), 
we obtain 



Pi 



<PL Pl 1 



h 



h 2 



dt'PL exp 



h 



{t-t'){l-P)L 



'l-P)L Px {t'). (A7) 



Using the definition for the Liouvillian superoperator L given above and assuming 

Tr env (6 n p env ) Tr env (c^p env ) 0, 



(A8) 



Eq. (A7) simplfies to 



h = - l h LsPS -J? J0 
where the self-energy superoperator kernel is 



dt'Zsit-t^psit'), 



(A9) 



E s (t) = Tr env [Ls-env exp(-iQLt/ h)Ls- env Pe 



(A10) 



and where Ls is the system part and Ls_ env the system-environment part of the full Liou- 
villian superator L. 

We now make the Born approximation, which amounts to dropping the interaction part 
Ls-env from the full L appearing in the exponential term of the kernel: 



£ s (t) « Tr cnv [L s _ eQv e~^ Ls+L -W h L s _ env p env ] 
= Tr \Lq e -i(L s +L cnv )t/hT i 

11 env \_ J— env° - LJ b—envycnv j > 



(All) 



where we can drop the Q projector as in the last line, a consequence of Eq. (A8). We 



assume that the system-environment interaction is sufficiently weak to justify making this 
Born approximation. 

Given the bilinear in operators form of the system-environment interaction Hamiltonian 
Hs-cnv [third and fourth terms in Eq. (26)] and using the following identity for any two 
operators A and B 

1 3 

Tt(AB) = - J2^(Aa v )Tr(Ba v ), (A12) 

u=0 
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where <7j, i = 1,2,3 are the Pauli matrices and <j$ is the identity matrix, the Born approxi- 



mation to the master equation ( A9 ) can be rewritten after some algebra as follows: 

ft 



Psit) 



l -[H s ,p s {t)]-±- 2 J dt'Tr cnv [BB(-f) Pcnv } |y, Y(—t')e~ iHst '^ h p s (t — t')e +iHst '^ h 



N 



h 2 

+h.c. 



(A13) 



In this expression, 0(t) = e lHst ' h Oe~ lHst ' h for some operator O, "h.c." denotes the hermitian 
conjugate of the preceding integral terms, and we have used the shorthand notations B = 
K n {bn + b n ) and = J2 n ^n\cn^ + + Cn ). Using the fact that the environment is in a 
thermal state to work out the environment correlation functions, we obtain: 



oo 

,.i, v ; dwJ(u)[coth(Phw/2)cos(ut)+ism(ut)] 
o 



Tr cnv [BB(t)p, 
where the bath spectral density is 

n 

and we have a similar equation for the correlation function involving C^: 

POO 

Tr env [C u) C u \t)p env ] = / dujJ {j \uj) [coth(/3/iu;/2) cos(c;t) + % sin(wt)] 

Jo 

with 

n 

For simplicity, we assume an "ohmic" spectral density with power law cutoff: 

A 2 

J(uj) = au 



(A14) 



(A15) 



w 2 + A 2 ' 



(A16) 



(A17) 



(A18) 



where A is the ultraviolet cutoff frequency. Substituting (A18) into (A14) and using contour 
integration to solve for the integral, we obtain: 

7raA 2 



Tr env [BB(t)p c 



e~ A * cot (J3KA/2) + -r— > - + ie~ M 

V 1 ' ph 2 ^ (2nn/f3hf -A 2 



(A19) 



with a similar expression for the correlation function (A16). 



From the form of the environment correlation function ( A19 ), we see that it decays rapidly 



to zero relative to the oscillator and TLS dynamical timescales, provided we assume that 
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the environment temperature satisfies fc^T ^> hu,E^\ where = \J (A^) 2 + (A ^ 2 



is 



the jth TLS energy level separation. Subject to this condition on the temperature, we can 



make a Markov approximation in master equation (A13) by setting t — > oo in the upper 
integration limit and expanding to first order in time the system's t'-dependent unitary 
evolution operator wherever it appears, i.e., e ±lHst ' ~ l±iHgt' '. The resulting Born-Markov 
master equation is 

Ps(t) = ~^[H s ,p s (t)} 

+§[Y 2 ,Ps(t)\ - ^[Y,{P Y , Ps (t)}] - *±[Y, [Y, Ps (t)\] + ^[Y, [Py,ps(t)\] 

1 N . N 



3=1 11 3=1 

1 N 

+^J2 A b )A 2WM j) ,Ps(t)}}, (A20) 

3=1 

where Py is the oscillator momentum, {•, •} denotes the anticommutator, and the R and / 
coefficients are the real and imaginary parts of the environment correlation function time 
integrals: 

POO 

R 1 = Re / dtTi cnv [BB(t) Pcnv ] , (A21) 
Jo 

POO 

R 2 = Re / dttTi cm [BB(t)p env ] , (A22) 
Jo 

POO 

= Re / dtTr env [C^C^\t)p cnv ] , (A23) 
Jo 

POO 

Rf = Re / dttTr env [C {j) C ij \t)p cnv ] (A24) 
Jo 



and analogously for the imaginary parts. In the Born-Markov master equation (A20), the 

11 term renormalizes the frequency of the oscillator, while the R 2 term is the so-called 
'anomalous diffusion' contribution. We shall neglect both terms, justified because of the 
assumed weak system-environment coupling and the above condition on the environment 
temperature. The remaining oscillator environment terms involving the I 2 and R\ cause 
damping and thermal diffusion of the oscillator, respectively. It is convenient to parametrize 

1 2 in terms of the energy damping rate 7 of the oscillator in the absence of the TLS: 
I 2 = Hm'y/2. The diffusion coefficient then becomes Ri = rwykBT, the expected form 
that follows from the fluctuation-dissipation theorem. The effect of the remaining three 
TLS environment terms are most straightforwardly understood by considering the coupled 
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moment equations of the three Pauli matrices in the absence of the oscillator. One finds 



for weak system-environment coupling that the R^ term can be neglected, while the R\ 
term causes damping/dephasing of the TLS and the diffusion term ensures that the 
moments decay to the thermal equilibrium state. Again, it is convenient to parametrize 
Ri in terms of the relaxation time T± of the jth TLS excited eigenstate in the absence of 



?0) 



the oscillator. From the moment equations, we obtain R^ = h?(E^ /A b ) 2 (4T] ) 1 and 
l!f = h/ (2k B T)R!f\ as follows from the fluctuation-dissipation theorem. In terms of these 



parametrizations, master equation ( A20J) becomes 



ps® = -l[H s ,Ps{t)}- l ^%{p Yl p s m- 1J ^^[YAy,psm 



N 

-E ! 



(j) 




^\ps{t)\\ 



{°i j) ,Ps(t)}\, 



(A25) 



— 8T? 

where we recognize in the first line the familiar master equation for a quantum Brownian 
oscillator in the large temperature limit ™ 



While the master equation (A25) is valid in the large temperature limit k B T Huj, 



it is desirable to investigate the system dynamics at low temperatures as well, such that 



kgT < fkv,A^\ In principle, a more involved analysis of Eq. (A13) with correlation re- 



lation expressions (A19) can yield a Markovian approximation that is valid at lower tem- 



peratures. However, a more direct way is simply to invoke the quantum version of the 
fluctuation-dissipation theorem, which for a Brownian oscillator amounts to making the 



following replacement in the diffusion term in (A25): 



tanh 



2k B r \2k B ry (A26) 

The resulting Born-Markov master equation describing the quantum Brownian motion of 
the oscillator alone is valid for temperatures k B T ^> 7ft. Given the assumed weak system- 
environment coupling, i.e., a large quality factor Q = ^> 1, oscillator dynamics can now 
be investigated at low temperatures such that k B T < hu. Analogously, we can make the 
replacement 



2k B T 



tanh 



gg) 
2k B T 



(A27) 



in the second TLS (diffusion) term of the master equation (A25). The resulting moments for 
the three Pauli matrices in the absence of the oscillator then decay to the correct quantum 
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thermal equilibrium state as required. We shall assume that making the replacements (A26) 



and (A27) in (A25) yield the Born-Markov master equation that is valid for temperatures 
ksT ^> jh,h/Tx, provided the interactions between the oscillator and TLS, as well as 
between the oscillator-TLS system and environment are weak, i.e., < hu,E^; Q > 1 
and E^Tt/h > 1. 
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